Accelerating cycle expansions by dynamical conjugacy 
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Periodic orbit tiieory provides two important functions — the dynamical zeta function and the 
spectral determinant for the calculation of dynamical averages in a nonlinear system. Their cycle 
expansions converge rapidly when the system is uniformly hyperbolic but greatly slowed down 
in the presence of non-hyperbolicity. We find that the slow convergence can be associated with 
singularities in the natural measure. A properly designed coordinate transformation may remove 
these singularities and results in a dynamically conjugate system where fast convergence is restored. 
The technique is successfully demonstrated on several examples of one-dimensional maps and some 
remaining challenges are discussed. 



I. INTRODUCTION 

Equilibrium statistical physics has been extremely 
successful, while accurate computation of physical 
averages in non-equilibrium systems remains a great 
challenge both theoretically and practically, due to 
the intrinsic difficulty of procuring the right statisti- 
cal weight of various states based on equations of mo- 
tion that govern system evolution [3| • From a dy- 
namical systems point of view, the statistical weight 
is proportional to the natural measure of states in 
the phase space, which often is non-smooth or even 
singular in a chaotic system [3 113 and thus flunks 
an accurate representation. Fortunately, periodic or- 
bit theory (POT) avoids a direct description of the 
possibly fractal measure by expressing phase space 
averages in terms of averages on periodic orbits or 
cycles and thus is a powerful way for reliable and 
accurate characterization 0, [l^ Ho, [13 of a nonlin- 
ear chaotic system. The associated cycle expansion 
of spectral determinants or dynamical zeta functions 
orders cycles in a hierarchical way such that dynam- 
ical averages are dominated by a few short cycles 
and the longer ones give decreasing corrections. For 
a uniformly hyperbolic system with finite symbolic 
dynamics, the corrections decrease exponentially or 
even super-exponentially P, [O, [HI- However, for 
non-hyperbolic systems, the convergence could be 
extremely poor, which severely limits the applica- 
tion of cycle expansions ■ 

Real physical systems are non-hyperbolic in dif- 
ferent cases. For example, in an intermittent sys- 
tem, typical trajectories alternate between regular 
and chaotic motions in an irregular way and thus 
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cause non-hyperbolicity [3, [ll|. A milder type of 
non-hyperbolicity is created by the strong contrac- 
tion at specific points such as homoclinic tangencies 
in the Hcnon map or critical points in 1-d maps 1, 2] . 
As a consesuence, the natural measure becomes 
non-analytic and thus the nice shadowing property 
among cycles that is necessary for fast convergence 
of cycle expansions fades out. Poles appear near the 
origin of the complex plane in the dynamical zeta 
function and the spectral determinant, which gives 
much trouble to a polynomial approximation as in 
the normal cycle expansion. To compute averages 
with fair accuracy, many cycles are needed, which 
usually requires unaffordable amount of resources. 
Thus, how to accelerate the convergence of cycle ex- 
pansions in the presence of non-hyperbolicity is a 
key problem in practice. 

Several accelerating schemes have been proposed 
based on the analyticity of the spectral functions. 
One idea is to identify and remove the poles that are 
near the origin and thus expand the radius of con- 
vergence. In P, H, 0] J the dominant terms in the tail 
of the expansion are estimated and summed up to 
approximately determine the leading pole. More ac- 
curate estimation is obtained by using Pade approx- 
imation, which is valid not only for compu ting the 
leading pole but also for seeking other ones |13|. IT8|. 
An interesting consequence of analyticity of the un- 
derlying dynamics is the existence of infinite sum 
rules which symptom strong correlations among pe- 
riodic orbits. These exact relations show signs of 
information redundancy embedded in the whole set 
of periodic orbits and can be utilized to accelerate 
the convergence of cycle expansions Tol. In certain 
cases, analyticity can be used to derive the spectral 
function with no resort to periodic orbits and thus 
implies a potential alternative route to the spectrum 
computation. But so far, it succeeded only for sev- 
eral very specific maps and hard to be generalized 



to other examples Q. 

In this paper, we employ a geometric picture of 
the cycle expansion to treat the convergence prob- 
lem for 1-d maps @, Q . Maps with critical points 
have a natural measure with singularities, due to 
the strong contraction around critical points. This 
contraction also deteriorates the dynamical shadow- 
ing between cycles of different lengths and thus leads 
to a slow convergence of cycle expansions. So, the 
singularity in the natural measure is an effective in- 
dicator of unbalance of cycle weights and signals a 
small radius of convergence. One idea for expedit- 
ing convergence is thus to identify and then clear out 
singularities in the natural measure. In the current 
paper, we achieve this by properly designing coordi- 
nate transformations such that the resulted conju- 
gate map has a natural measure with no singularity. 
The computation of dynamical averages in the orig- 
inal map can be efficiently done with counterparts 
in the conjugate map since the convergence is much 
improved in the new map. 

In the following, after a brief review of periodic or- 
bit theory in Section |lll we discuss in Section Hill the 
convergence of the dynamical zeta function and the 
spectral determinant for maps with critical points. 
A comparison between the two spectral functions is 
made and the importance of hyperbolicity for effi- 
cient calculation is emphasized. With a description 
of the geometric significance of the truncation in cy- 
cle expansions, our accelerating scheme is presented 
and tested on several examples. In Section IIVI we 
summarize the paper and discuss the existing prob- 
lems and possible directions for further investiga- 
tion. 



II. PERIODIC ORBIT THEORY 

More often than not, dynamical averages are con- 
veniently computed via time averaging, 

I N 

"■(xo) = —^^aix^), (1) 

i=0 

where xq — >■ xi — ^ • • ■ Xi — >■ ■■•xn is an itinerary 
generated by the map f{x) and iV is a large num- 
ber. Time averaging is easy to do but hard to 
achieve high accuracy. In the presence of non- 
hyberbolicity, its convergence is very slow and the 
result becomes unreliable. To better understand the 
dynamics for more efficient calculation, the phase 
space average (a) is introduced and has the nice 
property (a) = a(xo) in an ergodic system. Thus, 



a geometric picture is enabled to explore the averag- 
ing process, which will be explained in more detail. 

For an ergodic map /(x), when n — oo, under the 
map action any smooth initial measure will approach 
an invariant measure, called the natural measure. 
Formally, 

p{x)= hm / dy6{x-riy))poiy), (2) 

where p{x) is the natural measure and po{y) is an 
initial smooth measure. With the natural measure, 
the average (a) can be obtained by 

(a) = / a{x)p{x)dx , (3) 
Jm 

which is the most common way in statistical physics 
for computing averages. In a chaotic system, how- 
ever, the measure p{x), being often singular and sup- 
ported on a fractal set, is hard to obtain, which mo- 
tivates introduction of periodic orbit theory. 

For a map f : M ^ A4 and an observable a{x), 
we define the evolution operator £" o g[y) = 

C"^{y, x)g(x)dx for any function g(y). The kernel 

is 

C-{y,x)^5{y^ r{x))eP^\ (4) 

where j3 is an auxiliary variable and A" = 
X]fc=o ^{f^{x))- The average (a) or other dynamical 
properties can be conveniently obtained by virtue of 
the spectrum of C. Suppose the leading eigenvalue 
of C is e'*" , then we have 0] 

(.) = ^ . (5) 

Specifically, when /3 = 0, we have £"-{y,x) — 6{y — 
/"(a;)), which is the kernel of the so-called Perron- 
Frobenius operator. The escape rate of a dynamical 
system, which we denote by 7, can be obtained by 
computing the leading eigenvalue of this operator 

7 = -So- (6) 

The eigenvalues of the evolution operator £ can be 
detected with the help of the spectral determinant, 
which is related to the trace of £ and thus to the 
periodic orbits by the identity 

lndct(l - z/:) = Trln(l - z/:) . (7) 

For one-dimensional maps, detailed manipulation 
shows that [7] 

det(l-z/:)=exp(-5:5:---— ), (8) 
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where p denotes prime cycles which are not repeats 
of shorter ones. Ap is the stabihty eigenvalue of cycle 
p and Up is its length. In most classical computa- 
tions, we are only interested in the leading eigen- 
value. Obviously, the smallest positive zero of the 
above-defined spectral determinant is e~''° , the in- 
verse of the leading eigenvalue of £. In view of 
Eq. (O, we are able to compute dynamical averages 
with the spectral determinant. 

The leading eigenvalue can alternatively be ob- 
tained from a simpler spectral function — the dynam- 
ical zeta function 

l = [](l-ip),ip = -L,n.e/^-4''. (9) 



It can be proved that ^ is the zeroth-order approx- 
imation of the spectral determinant and has kept 
the smallest positive zero unchanged. Most often, 
however, they have different analytic properties. 

Practically, to evaluate zeros, we expand the spec- 
tral determinant or the dynamical zeta function in 
terms of power series and get polynomial approx- 
imation through truncation. The power series ex- 
pansion is one type of cycle expansion. For example, 
for the one-dimensional map with complete binary 
symbolic dynamics, the dynamical zeta function can 
be expanded as 

i = 1 - to - h - [{toi - hh)] 

— [(^001 " ^Ol^o) + (^011 — ^oi^i)] — • • • . 

(10) 

If we keep only the terms explicitly shown in 
Eq. pH)) . its cycle expansion is truncated at cycle 
length 3 and results in a polynomial in z of degree 
3. The linear term is the fundamental contribu- 
tion which gives the dominant part of the expansion. 
Higher order terms are curvature corrections which 
consist of contributions from prime cycles such as ^oi 
and from pseudo-cycles being combination of prime 
cycles such as toti. The cancelation between cycles 
and pseudo-cycles signals shadowing properties and 
smoothness of the underlying dynamics and results 
in an exponential decrease of the curvature correc- 
tions when uniform hyperbolicity is assumed. How- 
ever, there is no good cancelation in the presence 
of non-hyperbolicity and as a result cycle expansion 
converges very slowly. In the current paper, we are 
trying to restore this cancelation by dynamical con- 
jugacy under certain circumstances. 



III. IMPROVING THE CONVERGENCE OF 
CYCLE EXPANSIONS 

A. Notes on numerical computation 

To calculate dynamical averages, we need to build 
the truncated version of the dynamical zeta function 
and the spectral determinant. The detailed explana- 
tion for an efficient computation can be found in [7] , 
which is omitted here for brevity. With a truncation 
length TV, we drop out all the cycles longer than N. 
To study the convergence of cycle expansions, as an 
example, we will evaluate the escape rate and other 
dynamical averages with the truncated dynamical 
zeta function and spectral determinant. Also, we 
will check how the computational error of physical 
averages depends on the truncation length. How- 
ever, except for very few cases, we cannot obtain 
the exact average values. Therefore, we use aver- 
ages obtained with the truncation length Ny^nx + 1 
as the "exact" values when estimating errors with 
the truncation length no larger than A'max- Figures 
are plotted to show this dependence and the loga- 
rithmic scale is often used in the ordinate. 

The natural measure is computed by map itera- 
tions. We choose a random initial point and iterate 
it many times, usually 10^ if not specified otherwise. 
Then, by counting the times that the point enters 
a small interval, we get a probability distribution, 
which is a numerical approximation of the natural 
measure for an ergodic system. As the maps dis- 
cussed in this paper are all ergodic, though not very 
accurate in some cases, this method is simple for 
getting a rough picture of the natural measure. In 
addition, with Eq. ([T]), physical averages are easily 
computed with the iterations at the same time. 



B. Comparison of the dynamical zeta function 
and the spectral determinant 

The escape rate and other dynamical averages 
may be evaluated with the dynamical zeta func- 
tion or the spectral determinant. However, the con- 
vergence rates of the two methods are quite differ- 
ent, which is due to the difference in their radius 
of convergence in cycle expansions. To show this, 
we calculate the escape rate of the map f{x) = 
6x{l - x), M ^ [0,1], f : M M with the dy- 
namical zeta function and the spectral determinant. 
The results are listed in TABLE HI According to 
the table, both methods converge fast, thanks to 
the nearly perfect cancelation between prime and 
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0.831493012 


0.831492987487621 


6 


0.831492987 


0.831492987487621617307 


7 


0.8314929875 


0.8314929874876216173072762950 


8 
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0.8314929874876 




10 


0.83149298748762 





TABLE I. Escape rate obtained by the dynamical zeta 
function i and the spectral determinant det(l — zC) for 
the map fix) — 6x{l — x) on the interval [0, 1]. 
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(a)the error of the escape rate computed with the dynamical 
zeta function 



pseudo-cycles. Thus very accurate results can be 
obtained with only several short prime cycles. More- 
over, the results computed with the spectral deter- 
minant converge much more quickly than with the 
dynamical zeta function, implying a difference in 
their analyticity. 

FIG. [T] shows the dependence of the error in the 
computed escape rate on the truncation length N 
with the dynamical zeta function and the spec- 
tral determinant. It is clear that the logarithm of 
the error decreases linearly for the dynamical zeta 
function and super-linearly for the spectral determi- 
nant, which suggests an exponential and a super- 
exponential decrease in the error itself, respectively. 
The reason for this difference is that the spectral de- 
terminant is analytic over the whole complex plane, 
while the dynamical zeta function is only analytic in 
a region with a finite radius. Thus, the coefficient 
of the A^th-order term decreases super-exponentially 
with N in the spectral determinant, and exponen- 
tially in the dynamical zeta function. From a pure 
algebraic point of view, the way in which the coef- 
ficient of the A^th-order term decreases determines 
the convergence rate. 



C. The influence of hyperbolicity 

Having compared the convergence rate of the dy- 
namical zeta function and the spectral determinant 
for the map f{x) = Ax{l — x) with A = 6, we check 
how the value of A influences the convergence rate. 

First, we set A equal to 5 and repeat the above 
computation for the escape rate, the results are 
shown in FIG. [2j We see that though a little 
slower than the A = 6 case, both methods converge 
fast — the dynamical zeta function method converges 
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(b)the error of the escape rate computed with the spectral 
determinant 

FIG. 1. The error of the escape rate for the map f{x) = 
6x{l — x) computed with (a) the dynamical zeta function 
and (b) the spectral determinant. 



nearly exponentially and the spectral determinant 
exhibits a beautiful super-exponential convergence, 
just like what happened before. It looks as if the 
change of A had little effect on convergence. How- 
ever, if we set A equal to 4, a dramatic change hap- 
pens to the convergence. As shown in FIG. [3l the 
dynamical zeta function with cycles up to length 15 
gives an error of 10~^, while in the A = 5 case, the 
error is 10~^^. Moreover, the results obtained by the 
spectral determinant even lose the super-exponential 
convergence and exhibit only an exponential conver- 
gence. This phenomena implies that, in this special 
case, the spectral determinant may not be an entire 
function any more. 

If we increase A from 4 to 4.001, the convergence 
of the expansion improves dramatically, as depicted 
in Figure m When A — 4.001, the super-exponential 
convergence of the spectral determinant is restored. 
An apparent property that makes the map of A = 4 
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(a)the error of the escape rate computed with the dynamical 
zeta function 



(a)the error of the escape rate computed with the dynamical 
zeta function 
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(b)the error of the escape rate computed with the spectral 
determinant 



(b)the error of the escape rate computed with the spectral 
determinant 



FIG. 2. The error of the escape rate for the map f{x) = 
5x{l — x) computed with (a) the dynamical zeta function 
and (b) the spectral determinant. 



FIG. 3. The error of the escape rate for the map f{x) = 
41(1 — 2;) computed with (a) the dynamical zeta function 
and (b) the spectral determinant. 



different is that the height of the critical point faUs 
within the interval [0, 1], which is known to be the 
cause of the slow convergence j7|]. Here, we study 
this phenomena in great detail and will design a 
technique to counter its effect later. To show how 
the critical point influences the convergence of cy- 
cle expansions, we use the dynamical zeta function 
to calculate the escape rate for maps with a higher- 
order critical point. One general form of such maps 
is fix) = 1 - |2j; - 11'=^, a; e [0, 1] 

In FIG. [5l the profile of the map for different k 
is portrayed. As the order k increases, the top of 
the map profile becomes flatter and flatter. The er- 
ror of the escape rate obtained from the dynamical 
zeta function is displayed in FIG. [Hand FIG. [2Ia), 
where we can see that the convergence is poorer for 
maps with a higher-order critical point. To know 
why the dynamical zeta function and the spectral 
determinant flaw for maps with critical points, we 



must have a clear understanding of the nature of 
the approximations when we apply a truncation to 
the dynamical zeta function. 



D. The significance of the truncated dynamical 
zeta function 



As the number of all prime cycles is infinite, 
the truncation to the spectral functions is needed 
for an efficient computation. For example, for the 
unimodal maps discussed above, if we truncate at 
the shortest cycles, the dynamical zeta function is 
^ = 1 — to — ti. This truncated dynamical zeta func- 
tion leaves out all the curvature corrections and is a 
rough approximation for the original map. Now, one 
question can be asked: is there a linear map having 
^ — 1 — to ^ ti as its truncated first-order dynami- 
cal zeta function? A simple example is a piecewise 
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(a)the error of the escape rate computed with the dynamical 
zeta function 



FIG. 5. The graph of f{x) = 1 - (2x - 1)*, fc = 2, 4, 6. 
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(b)the error of the escape rate computed with the spectral 
determinant 



FIG. 4. The error of the escape rate for the map 
f{x) = 4.001a;(l — x) computed with (a) the dynami- 
cal zeta function and (b) the spectral determinant. 



linear map consisting of two branches, which is con- 
structed in this way: firstly, find the fixed points of 
the map along with the slopes at these points, then, 
draw line segments which pass the fixed points and 
are tangent to the graph of the map. If we want to 
construct a piecewise linear map which has a dynam- 
ical zeta function identical with that of the original 
map up to order A'^, we need the positions and slopes 
of all the periodic points of period not larger than 
A^. FIG. [7] displays such a piecewise linear map for 
the truncation up to length two, while the original 
map is /(x) = 4a;(l — x). 

For the piecewise linear map, the curvature cor- 
rections with order lager than TV are nearly zero due 
to its linearity, and therefore, we regard it as a proto- 
typed geometric model of the A^th-order truncated 
dynamical zeta function. The dynamical averages 
computed with the A'^th-ordcr truncated dynamical 
zeta function are very close to those given by the 



piecewise linear map. So, how well the piecewise lin- 
ear map approximates the original map determines 
the computational accuracy of the truncation. 

We know that the average obtained from the dy- 
namical zeta function is the phase space average 
(a) — Jj^ dxp{x)a(x), where p{x) is the natural mea- 
sure. For example, for the tent map, the natural 
measure is uniform everywhere. However, for the 
map with critical points, the natural measure has 
singularities as portrayed in FIG. [H] 

It's just these singularities that lead to the slow 
convergence in a cycle expansion, because the piece- 
wise linear map obtained from the truncated dy- 
namical zeta function can not capture well the nat- 
ural measure with singularities. Though the natural 
measure of the piecewise linear map gets closer and 
closer to that of the original map with increasing 
truncation length, it fails at the singularity. As a re- 
sult, the average obtained from the piecewise linear 
map, or equivalently, from the truncated dynamical 
zeta function for a non-hyperbolic system, does not 
converge as fast as in the uniformly hyperbolic case. 



E. Accelerating convergence in the presence of 
a critical point 

The singularity in the natural measure causes the 
slow convergence of the dynamical zeta function. A 
natural cure of this trouble is to clear out the sin- 
gularity by a coordinate transformation, which re- 
sults in a new map conjugate to the original one but 
without a singularity in the natural measure. Con- 
sequently, we are able to accelerate the convergence 
of the cycle expansion with the help of the conjugate 
map. 
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(a)the error of the escape rate for the map 
f(x) = l-{2x- 1)'^ 
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(b)the error of the escape rate for the map 
fix) = l-{2x- 1)6 

FIG. 6. The error of the escape rate for maps with a 
higher-order critical point computed with the dynamical 
zeta function. 




X 



FIG. 7. A piecewise linear map approximation of the 
map f{x) = 4a;(l — x). 



FIG. 8. The natural measure of the map f{x) = 4x{l — 
x). 



1. Clear out the singularities 

For maps with a critical point, such as f{x) = 1 — 
|2a; — l]*^, fc = 2,4,6,---, the probability distribution 
has an algebraic form in the neighborhood of the 
singularity, more explicitly 

P{^) ^ k-i near a; = , (11) 

where k is the order of the critical point. For ex- 
ample, the natural measure of the map f{x) — 
1 — (2a; — 1)^ has two singular points: and 1, with 
the probability distribution p{x) ~ near x = Q 

and pix) ^ ,} near x — 1. 

A coordinate transformation is needed which 
stretches the coordinate axis around the singular- 
ity in order to remove it. We are able to con- 
struct a homeomorphism h : M. ^ M. oi the form 
h{x) DC |a; — xq^^^ in the neighborhood of the sin- 
gularity X — xq to achieve this goal. For the map 
f{x) ^l-{2x- If = 4a;(l - x), [0, 1] ^ [0, 1], an 
appropriate transformation is h{x) = ;| arcsin -^/a;, 
which as you can see has the desired form in the 
neighborhood of and 1. 

With the coordinate transformation, the original 
map is changed to its conjugate. Denoting the origi- 
nal map by f{x) and the conjugate by g{x'), we have 
the relationship that / = o goh, or equivalently, 
g = ho f o /i-i. The map f{x) = 1 - (2x - 1)^ = 
4x(l — x) and its conjugate are displayed in FIG. O 

The conjugate map has no critical point any more, 
and the natural measure produced by the two maps 
are displayed in FIG. [TUl in which we can see that 
the natural measure of the conjugate map has no 
singularity any more. 
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(a)/(x) = 4x(l - x) 




(b)the conjugate map 

FIG. 9. (a)The map f{x) = 4x(l - x) and (b) its con- 
jugate map. 



2. The conjugate dynamical zeta function 



In an ergodic system, the dynamical average (a) 
can be obtained through the time averaging, as well 
as an average on the natural measure as discussed 
before. 

An ergodic map f : X ^ X and its conjugate 
g — h o f o h^^ : X' — > X' is related by the 
conjugacy h : X ^ X'. If we have the iteration 
f{xi) — Xi+i, under the conjugation h{xi) = x^, it 
becomes g{x'j) = a;-^^, which is to say, one trajectory 
{xi} for the map f{x) transforms into a trajectory 
{x'^} for the map g{x'). In particular, the cycle for 
the map f{x) is a cycle for the map g{x'). If the 
conjugacy h is piecewise smooth, a typical trajec- 
tory has identical weight in both coordinates, which 
suggests that the dynamical average can be similarly 
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(a)the natural measure of map f{x) = 4x(l — x) 
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(b)the natural measure of the conjugate map 

FIG. 10. The natural measure (a) of the map f{x) = 
4x(l — a;) and (b) of its conjugate map. 



computed with iterations of the map g{x'), 

N N 
i—i i—1 

(12) 

Eq. ([T2l) shows that computing the average (a) under 
the map f{x) is equivalent to computing the average 
(a o h~^) under the map g{x'). 

Based on the discussion above, it is natural to 
introduce a concept: the conjugate dynamical zeta 
function as follows. 

Suppose that the map f{x) and g{x') are conju- 
gate, with / — o g o h. The observable a for 
the map f{x) has a correspondent a o for the 
map g{x'). The dynamical zeta function ^ for f{x) 

and a, and ^ for g{x') and a o are said to be 

conjugate. We call h the conjugacy between ^ and 
j_ 

C'- ^ 

It's obvious that the dynamical averages obtained 
by Y and i are the same, since (a)/ = {a o h~^) g. 
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Typical forms for i and ^ are 



F. Several examples 



-c = np(i - tp). tp 



W (13) 
C - llpy^ ^ph - z |a;| ' 

where A'^ = EZi <^ih'\x',)) = = Ap. 

So, Ap is equal to A'^, no matter if h is diffeomor- 
phic or not. If h is diffeomorphic, the cycle stability 
eigenvalue Ap is also equal to Ap. Thus, we obtain 
i = ^. However, when h is not diffeomorphic, the 
cycle stability eigenvalue could be changed and so 
does the dynamical zeta function, i.e., ^ 7^ ^. This 
happens when we use a non-diffeomorphic coordi- 
nate transformation to clear out the singularities in 
the natural measure. 



The associated changes with the dynamical zeta 
function 



We mentioned in Section fill E II that to clear out 
the singularities, we need a coordinate transforma- 
tion h (X \x — Xq\^^'' near the singular point xq. The 
derivative of h has the form: 4^ oc r- So, the 

ax I 1 1 — -f ' 

dh/dx has a singularity aX x — x^. Thus, the conju- 
gacy h is not diffeomorphic at the singular point. If 
one periodic point happens to be singular, the stabil- 
ity eigenvalue of the periodic orbit will be changed. 

For example, for the map f{x) — \ — \2x— \\^ , the 
fixed point is a singular point of natural measure. 
So, the stability Aq for is changed to Aq , as shown 
below. 

In the neighborhood of the fixed point 0, the 
asymptotic form for / and h \s f ^ Kqx, x > Q 
and h^axk, where a > is a coefficient. We have 



g{x') = hofoh-\x') 



ho fixya'') - h{KQX^/a^) - A, 

(14) 

Hence stability Ao of for the conjugate map g is 
changed to Ay^ . 

So, in this situation, ^ 7^ ^- Nevertheless, for the 

map fix) = 1 — \2x — Ij'^, if the conjugacy is diffeo- 
morphic except at X = 0, 1, the only change in ^ is 

Ag = Ay^, as compared to ^. As the conjugate dy- 
namical zeta function ^ is computed for map g{x') 
whose natural measure has no singularity, the con- 
vergence for is much improved than for ^. Note 
that the exact functional form of the conjugacy is 
not essential as long as it has the right asymptotic 
form near the singular point. 



We proved that dynamical averages in a map can 
be computed with its conjugate map. For a map 
with critical points, we should find an appropriate 
coordinate transformation to clear out singularities 
in the natural measure. The conjugate map behaves 
much better in the sense that the singularity caused 
by critical points is eliminated and the convergence 
of the conjugate dynamical zeta function is accel- 
erated. Moreover, we do not have to know the ex- 
act functional form of the coordinate transforma- 
tion. What we do is change stability eigenvalues of 
specific cycles supported on the singularities of the 
natural measure and hence transform the dynamical 
zeta function to its conjugate. 

In the following, we apply our method to several 
examples. All these maps have critical points and 
therefore, produce a natural measure with singular- 
ity. 



1. The logistic map 

The logistic map f{x) — 4x(l — a;) ,x G [0, 1] has 
a critical point of order two. Its natural measure is 
singular at two points: a; = and a; = 1. Under the 
transformation h{x) = ^ arcsiny^, the singularities 
in its natural measure are cleared out, as shown in 
FIG.[TU1 Also, the conjugate map g{x') = ho f{x) o 
h^^ has no critical point any more, as depicted in 
FIG.Ii 

In the map g{x'), the stability eigenvalue of the 

fixed point is Ag — 2 where Aq = 4 is the stability 
eigenvalue of the corresponding point in the logistic 
map, with stability eigenvalues of all the other orbits 
remaining the same. 

The logistic map has an interesting property: the 
eigenvalue for any prime cycle except has an abso- 
lute value of 2", where n is the length of the cycle. 
However, the eigenvalue of is Ag = 4. After the 
coordinate transformation, the absolute value of the 
eigenvalue of any prime cycle of length n is 2", re- 
turning to the tent map case. 

Thanks to this interesting observation, the conju- 
gate dynamical zeta function under the conjugation 
h is just 



1 



(15) 



which is the same dynamical zeta function of the 
tent map. So, the escape rate of the logistic map is 
exactly 0. 
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Similar to the conjugate dynamical zeta function, 
we can write down the conjugate spectral determi- 

1/2 

nant, only with the change Ag = Aq' = 2. With 
the conjugate spectral determinant, the error of the 
escape rate decreases super-exponentially, as shown 
in FIG. [TT] However, if we apply the spectral de- 
terminant directly to the logistic map, the error of 
the escape rate decreases only exponentially with the 
truncation length TV, as shown in FIG.[3jb). So, by 
an appropriate coordinate transformation, we bring 
the super-exponential convergence back, which sig- 
nals that the influence of critical points wears out in 
the conjugate system. 

Certainly, the logistic map here is very special 
since it is exactly conjugate to a tent map [24| . To 
show the general applicability of the technique, in 
the following, we apply the method to several other 
maps for which no smooth conjugacy to the piece- 
wise linear map is known. 




FIG. 11. The error of the escape rate for the logistic 
map computed with the conjugate spectral determinant. 



2. The map f{x) = sin(7ra;) 

The map f{x) = sin(7ra;) , a; e [0, 1] has a critical 
point of order two. Similar to the logistic map, its 
natural measure has two singular points: x — Q and 
X = 1. The asymptotic form of its natural mea- 
sure near the singularity is p near a; = 

and p ^ ,} near x — 1. Under the coordinate 

transformation h[x) — ^ arcsin-^/a:, the singularities 
of the natural measure are removed, and we obtain 
the conjugate map g{x') = ho f{x) o hT^ . FIG. [12] 
portraits the map f{x) and g{x') where the critical 
point that exists in f{x) disappears in the map g{x'). 
FIG. [13] displays the natural measure produced by 



the map f{x) and g{x') respectively where the sin- 
gularity for f{x) at x = 0, 1 vanishes for g{x'). 




(a)the map f{x) = sin(7rx) 




(b)the conjugate map g(x') 

FIG. 12. The graph of (a) f{x) = sin(7ra;) and (b) its 
conjugate map g{x'). 

For the conjugate map g{x'), the stability eigen- 
value of the fixed point is changed from Aq = tt to 
Aq = TT^/^, while the eigenvalue of any other orbit 
doesn't change. We use the dynamical zeta func- 
tion and its conjugate to calculate the escape rate, 
(a;), (x^), {x^) of the map /(x) = sin(7rx). The re- 
sults are shown in TABLE [TT] with a cutoff of cycle 
length 10. Also included are the results obtained 
by direct time averaging, with lO'' iterations. From 
TABLE [n] we can see that the results obtained from 
the conjugate dynamical zeta function are by far the 
most accurate. 

FIG. [U displays errors in the averages obtained 
by the dynamical zeta function for f{x) = sin(7rx) 
and its conjugate with different truncation length. 
Although all the errors seem to decrease exponen- 
tially, the results from the conjugate dynamical zeta 
function decay much faster, indicating a great im- 
provement of the convergence. 
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(a)the natural measure of f{x) = sin(7rx) 



(a)the error of the escape rate 
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(b)the natural measure of the conjugate map g{x') 



(b)the error of (x) 



FIG. 13. The natural measure of (a) f{x) = sin(7ra;) and 
(b) its conjugate map g{x'). 





the dynamical 
zeta function 


the conjugate 
dynamical zeta function 


time 
averaging 


escape rate 


9 X 10"*' 


-3 X 10"" 




(x> 


0.47 


0.467962949 


0.468 




0.34 


0.34397492 


0.344 




0.28 


0.28394728 


0.284 



TABLE II. The escape rate, (x) , (x^) , (x^) for the map 
/ = sin(7rj;) computed with three different methods. 
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(c)the error of {x'^) 



3. The map f{x) = 1 - {2x - l)"* 

The map f{x) = 1 — (2x — 1)'* has a critical point 
of order four, with measure singularities at a; = 
and X — 1. The asymptotic form of the natural 
measure near the singularity is p ~ near a; = 

and p ' ^ 



■3- near x = 1 as shown in FIG. fUTa). 

(1-k)3 



To remove the singularities, an appropriate coordi- 
nate transformation is h{x) 



arccos(l — 2-^/1 — 

which has an asymptotic form h cx x^^^ near a; = 
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(d)the error of (x'') 



FIG. 14. The error of the escape rate, (x) , {x'^) , (x^) 
obtained by the dynamical zeta function(circles) for / = 
sin(7ra;) and its conjugate dynamical zeta function(stars). 



and h (X {1 — xY^^ near x = 1. Thus, we obtain the 
conjugate map g{x') = ho f{x) o h^^ . The map f{x) 
and its conjugate g{x') are depicted in FIG. [T5] Al- 
though f{x) has a very flat top, the peak of g{x') is 
acute. As a result, the natural measure of the map 
g{x') has no singularity, as exhibited in FIG. fW b). 




(a)the map f(x) = l-{2x- if 




(b)the conjugate map 

FIG. 15. The graph of (a) the map }{x) = 1 - (2x - 1)'' 
and (b) its conjugate map. 



The stability eigenvalue of of the conjugate map 
is Ag = Aq^^ — 8^/*. The convergence of the dy- 
namical zeta function for f{x) = 1 — {2x — 1)^ is 
even poorer than the logistic map. However, the 
conjugate dynamical zeta function continues to give 
a much accelerated convergence. The escape rate, 
{x) , (x^) , {x'^) obtained by the two different ways 



are shown in TABLE lllli with a truncation length 
10. It is worth mentioning that the direct time av- 
eraging becomes very unreliable in the current case. 
FIG. [T7] plots the errors in these averages with differ- 
ent truncation length. We can see that the conjugate 
dynamical zeta function converges much faster. 
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(a)the natural measure of f{x) = 1 — (2x — 1) 
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(b)the natural measure of the conjugate map 

FIG. 16. The natural measure of (a) f{x) = 1 — (22: — 1)** 
and (b) its conjugate map. 





the dynamical 
zeta function 


the conjugate 
dynamical zeta function 


escape rate 


2 X 10"^ 


-2 X 10"^ 


{x) 


0.45 


0.4475860 




0.36 


0.3601271 


{x^} 


0.32 


0.31801265 



TABLE III. The escape rate, (x) , (x^) , (a;^) for the map 
f{x) = 1 — {2x — 1)^ computed with two different meth- 
ods. 



I The map f{x) = 1 - (2s - 1)" 

The map f{x) = 1 — {2x — 1)^ has a critical point 
of order six, and therefore causes an even worse 
convergence for the dynamical zeta function. The 
coordinate transformation we use is h{x) = 1 — 

arccos(l — 2(1 — a: 3 ) ) 1.1 • - / /\ 

— — and the conjugate map g[x j — 

h o f{x) o h^^ has no critical point any more. Thus, 
the singularities of natural measure are removed. 
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(a)the error of the escape rate 
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(b)the error of (x) 
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FIG. [T^ shows the graph of the map /(x) and g{x'). 
The natural measure of the map f{x) and g{x') is 
depicted in FIG. [121 obtained with 10* iterations 
in this case. We see that both numerical measures 
fluctuate, which implies that time averages based on 
iterations can't reach a high accuracy for this map. 
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(a)the map f{x) = 1 - {2x - 1)^ 




(b)the conjugate map 

FIG. 18. The graph of (a) the map f{x) = 1 ~ {2x ~ if 
and (b) its conjugate map g{x'). 



(c)the error of (x^) 
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(d)the error of (x^) 



For the conjugate dynamical zeta function, the 
only difference from the original one is that, the sta- 

— - 1 

bility eigenvalue of is changed to Aq = Aq — 12<i . 

The values of the averages obtained by the two dy- 
namical zeta functions, with a truncation length 10, 
are listed in TABLE IIVI while the dependence of 
the computational errors on the truncation length is 
portrayed in FIG.[2D1 We can see that the conjugate 
dynamical zeta function converges much faster than 
the original zeta function, which provides evidence 
that clearing out the singularity in natural measure 
helps accelerate the convergence. 



FIG. 17. The error of the escape rate, (x) , (x^) , {x"^) 
obtained by the dynamical zeta function(circles) for 
f{x) = 1 — {2x — l)** and its conjugate dynamical zeta 
function(stars). 
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(a)the natural measure of the map f{x) = 1 — {2x — 1) 



(a)the error of the escape rate 
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(b)the natural measure of the conjugate map g{x') 



(b)the error of (x) 



FIG. 19. The natural measure of the map f{x) 
(2x — 1)" and its conjugate map g(x'). 





the dynamical 
zeta function 


the conjugate 
dynamical zeta function 


escape rate 


5 X 10"^ 


-2 X lO"'' 


{x) 


0.4 


0.40232 




0.34 


0.332921 




0.31 


0.30027 



TABLE IV. The escape rate, (a;) , (x^) , {x^} for the map 
f{x) = 1 — {2x — 1)^ computed with two different meth- 
ods. 



5. A map with three measure singularities 

The maps we have studied so far all have two mea- 
sure singularities: x = and x = 1. In this section, 
we turn to a map with three measure singularities, 
the graph of which is shown in FIG.I^Tla). The exact 
functional form of the map is f{x) = sin(^(l — x)), 
where a = 1.3156445888.... It has a critical point of 
order two, and has a nice property: /(O) — Xf, where 
Xf is the unique fixed point of the map. The natu- 
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(c)the error of {x'^) 
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(d)the error of (x'') 



'^^ FIG. 20. The error of the escape rate, (x) , {x'^) , {x^) 
obtained by the dynamical zeta function(circles) for 



fix) = l-{2x 
function(stars). 



1) and its conjugate dynamical zeta 



ral measure of / has three singularities: x = 0, a;/, 1. 
The asymptotic form of the natural measure near 
the singularity is: p ^ -k= near x = 0, p ~ , ^ 

near x — Xf and p ^ ^ near x = 1, as shown 

in FIG. [22r a'). To clear out the singularities, we 
use the coordinate transformation h{x) as depicted 
in FIG. [231 which stretches the coordinate around 
the singularities. The conjugate map g{x') is nearly 
a piecewise linear map, as shown in FIG. [2lT b). 
The natural measure of the map g{x') is depicted 
in FIG.mjb), which has no singularity any more. 




(a)the map with three measure singularities 




(b)the conjugate map 

FIG. 21. The graph of (a) the map with three measure 
singularities and (b) its conjugate map. 

For the conjugate dynamical zeta function, the 
stability oixj should be changed to \ Axf \' — lA^fl"^ ■ 
Again, we use the original and the conjugate dynam- 
ical zeta function to calculate averages. The results 
are listed in TABLE |Vl with a cutoff of cycle length 
20. The errors in the computation are plotted in 
FIG.[24l Note that the binary symbolic dynamics is 
not complete in the current example. The number 
of cycles get much reduced compared with the full 
symbolic dynamics case. 
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(b)the natural measure of the conjugate map 

FIG. 22. The natural measure of (a) the map with three 
measure singularities and (b) its conjugate map. 



By clearing out the singularities in the natural 
measure, the convergence is accelerated a lot. So, in 
this case, the conjugate dynamical zeta function is 
still an effective way to acquire averages with high 
accuracy. 





the dynamical 
zeta function 


the conjutate 
dynamical zeta function 


escape rate 


5 X 10"* 


-1 X 10"^° 


(x) 


0.601 


0.601895610 




0.453 


0.453165976 




0.36 


0.364669939 



TABLE V. The averages for the map with three measure 
singularities computed with two different methods. 
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FIG. 23. The conjugacy h{x) for the map with three 
measure singularities. 



A map with measure singularities on a period-2 
orbit 



The conjugation method may be applied to sys- 
tems with measure singularity on longer orbits. As 
an example, we study a map with measure singular- 
ities on a period-2 orbit. The functional form of the 
map is still f{x) = sin(-(l — x)), but with a dif- 
ferent value a = 1.10263451544766.... FIG. [^a) 
portrays the graph of the map. For this map, 
/(O) = Xa, /^(O) = Xb, where Xa and xi, make a 
period-2 orbit. The natural measure of the map 
is shown in FIG. I^STa) and has four singularities: 
x = 0, Xa, ccf,, 1. An appropriate conjugacy h(x) 
is plotted in FIG. [571 which stretches the coordi- 
nate around the four singularities. The conjugate 
map g(x') and its natural measure are plotted in 
FIG.Iinib) and FIG.[26lb) respectively. For the map 
g{x'), the singularities has been removed — the situ- 
ation that we have experienced many times. 

We change the eigenvalue of the period-2 orbit and 
obtain the conjugate dynamical zeta function — ^just 
as what we have done. The conjugate dynamical 
zeta function gives interesting results in the calcu- 
lation of the escape rate. FIG. [28] shows the error 
in the escape rate computed with the original and 
the conjugate dynamical zeta function. In general, 
the results obtained with the conjugate dynamical 
zeta function converge exponentially and uniformly, 
faster than with the original dynamical zeta func- 
tion. However, it's not totally true here. The conju- 
gate dynamical zeta function doesn't accelerate the 
convergence as effectively as before. So, why the con- 
vergence for the conjugate dynamical zeta function 
is not so good even if we have removed the measure 
singularities? 



(a)the error of escape rate 
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FIG. 24. The error of the escape rate, {x) , (a;^) , (x^) 
obtained by the dynamical zeta function(circles) for the 
map with three measure singularities and its conjugate 
dynamical zeta function(stars) . 



(a)the map with measure singularities on a period-2 orbit 



(a)the natural measure of map / 
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(b)the conjugate map 



(b)the natural measure of the conjugate map 



FIG. 25. The graph of (a) the map with measure singu- 
larities on a period-2 orbit and (b) its conjugate map. 



FIG. 26. The natural measure of (a) the map with mea- 
sure singularities on a period-2 orbit and (b) its conju- 
gate map. 



If we stare at FIG. [25t b). the graph of the con- 
jugate map, we find one special point at which the 
slope is infinite. It's just this point that causes the 
slow convergence. We regard this infinite slope map 
as "super-hyperbolic" . The super-hyperbolic map 
can slow down the convergence. To illustrate this 
point, we use the dynamical zeta function to calcu- 
late the escape rate of a super-hyperbolic map f{x), 



f{x) = { 



2x xe [0, 1/2] 

V2-2x xe [1/2,1], 



(16) 



where /'(I) = oo. We compare the results obtained 
from this map and the logistic map in FIG. 
Both results are computed with the original dynam- 
ical zeta function. We can see that the convergence 
rates are similar for these two maps. So, the super- 
hyperbolicity is harmful to the convergence, just like 
the non-hyperbolicity. To understand this, we recall 
the fact that the average obtained by the truncated 
dynamical zeta function is nearly identical to the 
one computed with the corresponding piecewise lin- 



ear map. So, if the piecewise linear map can ap- 
proximate the original map very well, the average 
obtained would be quite accurate. However, for the 
super-hyperbolic map, there exists a point with in- 
finite slope, which means that the value of the map 
changes extremely unevenly near that point. So, a 
lot of cycle points are needed near that point to get 
a fair approximation just like near the critical of a 
non-hyperbolic map. Thus, the slow convergence of 
the super-hyperbolic map is expected. 

Based on the discussion above, we can see that the 
singularity in the natural measure is not the only fac- 
tor which influences the convergence of the dynami- 
cal zeta function. For the map with measure singu- 
larities on a period-2 orbit, clearing out the measure 
singularities doesn't improve the convergence much. 
In fact, the direct and essential factor to determine 
the convergence of the dynamical zeta function is the 
cancelation between prime cycles and pseudo-cycles. 
For the super-hyperbolic case, the cancelation is not 
good even if the measure singularities do not exist. 
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FIG. 27. The conjugacy h{x) for the map with measure 
singularities on a period-2 orbit. 



FIG. 29. The error of the escape rate for the super- 
hyperbolic map and the logistic map. 
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FIG. 28. The error of the escape rate by the orig- 
inal(circles) and the conjugate(stars) dynamical zeta 
function for the map with measure singularities on a 
period-2 orbit. 



So, in the super-hyperbolic case, how to incur fur- 
ther cancelation remains a challenging problem. 



poles of the spectral function by classifying singu- 
larities in the natural measure and removes them 
through a coordinate transformation. 

We test our method on several maps, i.e., f{x) = 
1 - \2x - k = 2,4,6, /(x) = sin(7rx) which 
have only one critical point and complete binary dy- 
namics, and the map which has three measure sin- 
gularities. For these maps, the conjugate dynam- 
ical zeta function converges much faster than the 
original zeta function. Also, the conjugate spectral 
determinant restores the super-exponential conver- 
gence. However, when we treat the map with mea- 
sure singularities on a period-two orbit, we find that 
the conjugate dynamical zeta function does not con- 
verge as fast as expected. Analysis shows that the 
super-hyperbolicity of the conjugate map leads to 
this slowing-down. Further study is needed to elim- 
inate this nuisance. 



IV. CONCLUSION 

The central idea of this paper is that by clear- 
ing out the singularities in the natural measure we 
may accelerate the convergence of cycle expansions. 
Maps with critical points produce natural measures 
with singularities and show bad convergence in the 
expansion calculation. With appropriate coordi- 
nate transformation, the resulted conjugate map 
produces no singularity in its natural measure. To 
calculate dynamical averages, we use the conjugate 
spectral function, the cycle expansion of which is 
greatly accelerated due to the removal of the sin- 
gularities. Essentially, the method locates leading 



In this paper, we use one-dimensional maps as 
examples to demonstrate our accelerating scheme. 
How to generalize it to higher dimensions or to flows 
requires further investigation. Even in the 1-d case, 
we can can only treat maps with symbolic dynam- 
ics being subshifts of finite type. If the genealogy 
sequence of the critical point is not essentially pe- 
riodic, there may exist a natural boundary in the 
complex plane for the dynamical zeta function on 
which singular points are dense. In this case, the 
natural measure is singular on a dense and count- 
able set [l^l. It seems not possible to expand the 
radius of expansion by analytic continuation. Novel 
techniques have to be invented to achieve accelerated 
convergence in this case. 
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